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ABSTRACT 

We obtained global transonic solutions representing optically thin advection- 
■ dominated accretion flow by solving the full set of differential equations describing 

such systems. We found that far from the sonic point self-similar solutions are an 
excellent approximation of the global flow structure if the accretion rate is well below 
^ . the maximum value above which no optically thin solutions exist. 
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Accretion on to a compact object powers many astrophysical systems such as active galactic 
nuclei, X-ray binaries, and cataclysmic variables. It contributes to the energy sources of young 
stars (e.g., Bertout, Basri & Cabrit 1991) and plays a role in the physics of the nuclei of 'ordinary' 
galaxies (Fabian & Rees 1995) and of our Galactic Center (Narayan, Yi, Sz Mahadavan 1995). 
Until recently it was customary to assume, that the accretion power is radiated away with an 
efficiency of ~ GM/(rc 2 ) (where G is the gravitational constant, c is the speed of light, M is the 
mass of the central object, and r the radial distance from it), which is ~ 0.1 close to the surface 
of a neutron star or a black hole. In other words it was assumed that most of the gravitational 
energy released through viscous dissipation say, was locally radiated away from the accretion disk. 
This condition is very well satisfied for a geometrically thin accretion disk (e.g., Frank et al. 1992). 

Accretion flows in which radiative efficiency is very low and most of the heat gets advected 
into the central object have been studied by Begelman (1978, 1979) for spherical flows and by 
Begelman & Meier (1982), Abramowicz et al. (1986; 1988) for accretion disks. These disks were 
no longer geometrically thin (they were called 'slim disks' by Abramowicz et al. 1988) but were 
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optically thick and corresponded to very high, super-Eddington accretion rates. The application 
of these models to real objects have not yet been found. 

The interest in advection-dominated flows has been recently revived by the work of Narayan 
k Yi (1994, 1995a, 1995b), Abramowicz et al. (1995), Chen et al. (1995) and Chen (1995) 
in which optically thin solutions were found and studied. However, the main reason for the 
renewed interest in advection-dominated flows (e.g., Fabian k Rees 1995; Mineshige 1996) was the 
successful application of the model to the properties of the Galactic Center (Narayan et al. 1995) 
and the quiescent soft X-ray transients (Narayan, McClintock, k Yi 1996; see also Lasota 1996, 
Lasota, Narayan, k Yi 1996b). Recently Lasota et al. (1996a) applied advection-dominated flow 
models to the under- luminous active galactic nucleus NGC 4258. 

Despite of their successful applications, models of optically thin accretion flows need to be 
improved, at least in order to test the validity of the simplifying assumptions that are often made. 
First, the description of the microphysics involved requires improvements and completion. For 
example, pair creation has been neglected until now, but it has been recently included into the 
model (Bjornsson et al. 1996; Kusunose k Mineshige 1996). Second, the radial flow structure 
has been described in a simplified way. On the one hand Abramowicz et al. (1995) and Chen 
(1995), assumed that the flow is Keplerian (i.e. radially subsonic), but on the other hand, Narayan 
and collaborators use a self-similar solution which is significantly sub-Keplerian and subsonic 
everywhere. As shown by Chen et al. (1995) both type of models give the same general structure 
of the solutions but they differ in details. However, in the case of accretion on to a black hole 
the flow must be transonic (e.g., Abramowicz et al. 1988) so that neither of the descriptions is 
globally correct. 

The self-similar solution may provide an excellent description of an advection-dominated flow 
far from the flow boundaries as it was shown by Narayan k Yi (1994) in the case of a settling 
subsonic solution (i.e. accretion on to a non-magnetized star). However, the structure of the flow 
in the inner and outer regions may have an important influence on the general properties of the 
modeled objects: high energy radiation is emitted primarily close to the black hole whereas optical 
and UV radiation will originate mostly in the outer regions where the flow may form a Keplerian 
disk. 

In this paper we solve the problem of a transonic, advection-dominated flow on to a black hole 
and compare the solutions obtained with the local Keplerian and the self-similar solutions. The 
paper is organized as follows. In Section 2 we present the formulation of the problem. In Section 
3 we obtain a self-similar solution of the Narayan k Yi (1994) type which will be compared with 
the general transonic solution. Section 4 contains the results of the numerical calculations. Finally 
in Section 5 we discuss the differences between our solutions and the self-similar and Keplerian 
approximations. Section 6 contains the conclusions. 
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2. FORMULATION 



We consider an accretion disk which is axisymmetric, non-self-gravitating, optically thin and 
geometrically slim (H/r ^ 1) so that it can be described by vertically integrated equations. Here r 
and H are respectively the disk radius and the half-thickness. We will assume that the Keplerian 
rotational angular velocity is Qk = \J GMjrir — r g ) 2 implied by the pseudo-Newtonian potential 

= —GM/(r — r g ) (Paczyhski & Wiita 1980), where M is the mass of the central object and 
r g = 2GM/c 2 is the Schwarzschild radius. Equations describing a stationary accretion flow are 
summarized below. Specifically, the continuity equation is 

M = -2tt7tv£, (2.1) 

where M, v r , and £ = 2Hp are the mass accretion rate, the radial velocity and the surface density 
at a cylindrical radius respectively. Here p is the mid-plane density of the disk. The equation of 
motion in the radial direction is 

-^ + ^-^)r-~=0, (2.2) 

where Q, is the angular velocity and P is the vertically integrated pressure written, in a form 
compatible with the optically thin assumption, as 

V"7?T 

P = 2Hp= — , (2.3) 

where T and p are the mid-plane temperature and pressure respectively, ji is the mean molecular 
weight, which for cosmic abundances is equal to 0.62, and 1Z is the gas constant. The angular 
momentum equation can be integrated once with respect to the radius to yield the relation 

r 2 ^ = v r (£-£ m ), (2.4) 
dr 

where v is the kinematic viscosity and £ — £lv 2 is the specific angular momentum. Here £i n is 
the specific angular momentum lost from the disk, so that the viscous torque vanishes at the 
inner edge of the disk. The energy conservation equation is expressed by the balance between the 
local viscous heating, Q + , the local radiative cooling, Q_, and the global heat transport (radial 
advection), Q a dv It is expressed as 

Q+ = Q-+ Qadv (2.5) 
The viscous heating rate per unit area is given by 

2 

dr J 

The advection cooling rate is taken in a form (e.g., Chen k, Taam 1993): 



Q + (r^Y. (2.G) 



n - ^ VrP 



4-3/3dlnT o/fldjng 
r 3 -ldlnr 1 PJ dlnr 
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where (3 = P g /P is the ratio of the gas pressure to the total pressure, ^ = 
1 + (4 - 3/3) (7 - l)/[/3 + 12(7 - l)(/3 - 1)], and 7 is the ratio of specific heats. We will 
take (3 = 1 and 1^3 = 7. One should note that, due to different approaches to the vertical 
integration, the resulted Q a dv is different from the one of Narayan &: Yi (1994) as explained in 
Chen et al. (1995). The dimensionless advection factor £ characterises the entropy gradient. In 
our case, for a self-similar solution £ = 1 (see more discussion later). 

We shall assume that the local radiative cooling is provided by optically thin thermal 
bremstrahlung with emissivity (erg s _1 cm" 2 ), 

Q- = Qbrem = 1-24 X lO^Hp 2 ^ 2 . (2.8) 

Finally, we use a standard a- model viscosity prescription (Shakura & Sunyaev 1973): 

v = ^ac s H, (2.9) 

where a is a constant, c s = y/P/T, is the local sound speed, and H = c s /£Ik is the local scale 
height which is also the half-thickness of the disk. 

The above equations can be combined to reveal the presence of a sonic point. In particular, 
the equation of motion in the radial direction can be cast in the form 



r dty N 
v r dr D ' 



(2.10) 



where 



N = (tf K - tf)r 2 - Cl + (Q+ - Q-) , (2.11) 

ZjV r 



and 

D = C 2 S - v 2 r . (2.12) 
Here C s is the adiabatic sound speed (c s is the isothermal sound speed) defined as 

C 2 S = 7 f ■ (2-13) 

The Mach number is defined as: 

M = \v r \/C s . (2.14) 

The vanishing of both N and D at the sonic point provides the regularity conditions required 
for a transonic solution of the flow structure. Furthermore, the differential equation for the 
temperature can be expressed as 

lg- (1 . T) (" L^-V (2.15) 
T dr \ D Pv r V ; 
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In the steady state accretion disk equations, the specific angular momentum, £j n , is an 
eigenvalue of the problem which is adjusted in such a way that a transonic solution is found at 
the sonic point. The equations are solved numerically using a relaxation technique subject to the 
regularity condition at the sonic point (Chen & Taam 1993). The initial trial solution is obtained 
from the integration of the equations (2.10) and (2.15) by assuming dlnO/dlnr = dlnf^/dlnr 
(see Paczyhski & Bisnovatyi-Kogan 1981). The initial starting point for the integration is located 
away from the outer boundary used in the relaxation procedure. The exact location of the starting 
point does not affect the solution inside the outer boundary as long as it is sufficiently far away. 
Two kinds of advection-dominated solutions are used at the starting point. One is the Keplerian, 
advection-dominated solution of Abramowicz et al. (1995) and the other one is the self-similar 
solution (Narayan & Yi 1994). The trial solution of v r and T at the outer boundary is then used 
as the boundary value in the relaxation procedure (with dlnf£/dlnr ^ din fix/din r). Typically, 
the method allows to converge to a solution after a few iterations. 

At the sonic point, v r and T are calculated by requiring that both the numerator and the 
denominator in equation (2.10) vanish. Their derivatives are calculated by applying L'Hospital's 
rule. Defining y = din |v r |/dlnr at the sonic point, a quadratic equation for y can be derived using 
equation (2.10) and (2.15). For example, in a polytropic approximation, it reads: 

(7+l)y 2 + 2( 7 -l)y + 2a + 7 -l = 0, (2.16) 

where, 

„_ r d(^-^) (217) 



2C| dr 

(Note the difference with Chen & Taam (1993) where the last term 7 — 1 is missing). If y has no 
real solution, the sonic point is classified as a spiral point, which is not physical and is discarded. 
If the two solutions have different sign, the point is of saddle type, and if they have the same sign, 
the point is of nodal type (see Chen & Taam 1993 and references therein). 

For definiteness, we assume that the black hole has a mass of 1OM and the inner numerical 
boundary of the disk is chosen to be 1.5 r g . The outer numerical boundary of the disk is chosen 
to be typically about 1000 r g , but larger in models in which we discuss the effects of boundary 
conditions. The results do not depend on the spatial resolution which is typically Ar/r ~ 0.01. 



3. ANALYTIC SOLUTIONS 



In some cases the differential accretion disk equations can be simplified to an algebraic system. 
For example, in the limit of Keplerian disk models (where all the pressure-gradient and inertial 
terms are neglected), one obtains the Shakura & Sunyaev equations and the corresponding analytic 
solutions. Narayan & Yi (1994) obtained self-similar solutions describing advection-dominated 
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optically thin accretion flows. For the convenience of describing the numerical results in the next 
section, we derive the self-similar solution with the formulae outlined in §2. It should be mentioned 
that since the pseudo-Newtonian potential is used here, a self-similar solution is adequate only for 
large radii where the Keplerian rotation law can be approximated as fix = \/GM/r 3 . First we 
assume that each physical quantity X can be expressed as (Narayan & Yi 1994), 

X = Cr N , (3.1) 

where C = C(M,a). The power index N, which may be different for each physical quantity, is 
constant. For the advection-dominated optically thin accretion flows considered here, the power 
indices for v r , f2, and T are —1/2, —3/2, and —1/2 respectively. The solutions are: 

V r = ~ 5 t 2g 9^KT, (3.2) 

^t/^W (3.3) 
g(n K r) 2 , (3.4) 




M = ^/(5/2 + e')fl/7, (3.5) 

where, 

fe = e = — , (3.6) 

7-1 



«-V 1 + (6T5T*- L (3 ' 7) 

Here, / is the fraction of the advective cooling, and we have assumed / = 1 for the self-similar 
solutions. Note that e is different from that of Narayan & Yi (1994). This is due to a different 
definition of the advection term that was mentioned above. Note that, in the self-similar solutions, 
the dimensionless advection factor £ is determined by 7 only: £ = 1/(7 — 1) — 1/2. Therefore, 
£ = 1 for 7 = 5/3 independent on the value of / which one may have assumed. 

All other quantities such as p and H can then be calculated accordingly. The self-similar 
solution has its limitations since no boundary condition has been taken into account. It has a 
constant Mach number and therefore the transonic region can not be described. However, as it 
will be shown later, it describes correctly the true solution asymptotically at large radii if the flow 
is advection-dominated. 



4. GLOBAL NUMERICAL SOLUTIONS 



- 7- 



We have solved the full set of accretion disk differential equations numerically, using the 
method described in the previous section. The calculated models are summarized in Table 1 
in terms of the two dimensionless input parameters: mass accretion rate rh = M/Mg, and the 
viscosity parameter a. Here Me = ^.ixGM / '(cn es ) is defined as the Eddington accretion rate with 
K es = 0.34. Table 1 gives the eigenvalue £i n [in the unit of ^(3r g )] of the solution, the location of 
the sonic point r s , and the topological type of the sonic point. 

The form of the solution depends on the values of a and rh. Abramowicz et al. (1995), 
Narayan & Yi (1995ab) and Chen et al. (1995) found that for viscosity parameters a < a cr \t, 
where a cr i t is a critical value depending on r and on the assumed microphysics (see Chen et al. 
1995, Bjornsson et al. 1996), there exist a maximum accretion rate above which no optically thin 
solution exists. In the case of bremsstrahlung radiative cooling one gets (Abramowicz et al. 1995): 

m max » 1.7 x lOV/r^VVr 2 (4.1) 

For very small a, since m max oc a 2 r -1 / 2 , the optically thin solutions cease to exist for some 
r > r out , unless rh is very low. For this reason we could find global solutions only for not very 
small a or for very low accretion rates. 

We show in Figure 1 the results for a fixed value of a = 0.1 and two accretion rates: 
rh = 10~ 5 (Model 1, heavy solid line) and rh = 10 -2 (Model 2, heavy dashed line). The values 
of accretion rates are well below m max so that both solutions represent typical optically thin 
advection-dominated accretion flows (ADAFs). The form of the specific angular momentum is 
similar to that of the optically thick transonic accretion disks (see Abramowicz et al. 1988; Chen 
&; Taam 1993). The solution is everywhere sub-Keplerian because of the high value of a (see 
Abramowicz et al. 1988). The sub-Keplerian character of the flow is also seen in the absence of a 
pressure maximum (note that the quantity plotted is the pressure p, not the vertically integrated 
pressure P which always reach a maximum point independent on a). £ is everywhere close to 1 
and the self-similar solution is a good approximation for r > 30 r g . The rh = 10~ 2 solution is not 
as advectively dominated as the rh = 10~ 5 one but it is still a robust ADAF. £ is practically the 
same for both accretion rates confirming the analysis of Chen (1995). 

The dependence of the properties of the solution on a is illustrated in Figure 2. Here, the 
accretion rate is fixed as rh = 10" 5 and we consider two values of the viscosity parameter, a = 0.01 
(Model 3) and a = 0.001 (Model 4). For such low values of a the specific angular momentum 
becomes super-Keplerian in the inner disk regions. Accordingly a pressure maximum is present. 
For the higher value of a the flow is advection-dominated and £ differs from 1 only in the innermost 
parts of the flow. The character of the inner flow for both a = 0.01 and a = 0.001 is different 
from that of a = 0.1. In the first case (small a), the angular momentum i in the transonic region 
becomes almost constant with radius so that the viscous torque is very weak. This is therefore a 
'relativistic Roche-lobe overflow' type, pressure driven, accretion (see Abramowicz et al. 1988 for 
a detailed explanation). In the second case, there is a gradient in t and so a viscosity is needed 
for the mass to be accreted. This type of accretion could be therefore called 'viscosity dominated'. 
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One should also note that the sonic radius moves inwards with decreasing a. 

For a = 0.001 the situation is still different. Only the form of the angular momentum as 
a function of r is similar. The solution however ceases to be advection-dominated at r « 30 r g 
and £ is nowhere well represented by £ = 1. This is because, for this small a, the accretion rate 
m = 10~ 5 is near the maximum value for these radii, so the local cooling becomes important. The 
self-similar approximation is still a good representation of p for r > 30 r g but since £ is nowhere 
constant a self-similar solution is not a globally valid representation of the flow. It is seen that the 
local cooling tends to reduce £. In fact, for a radiative cooling dominated optically thin Keplerian 
disk, £ is negative (Wandel & Liang 1991). 

Finally, we present a case (see Fig. 3, Model 5) in which the sonic point is of nodal type 
(in all previous solutions the sonic point was of saddle type). In this case a = 1 and m = 10~ 3 . 
Because of the large value of a the flow is very sub-Kcplcrian and the sonic point is far away from 
the black-hole. £ is always bigger than 1 but all the particularity of this solution are due to the 
high viscosity, the type of the sonic point has no obvious manifestations. 



5. DISCUSSION 



In the inner transonic region the optically thin ADAFs have properties similar to those of 
the optically thick, radiation-pressure dominated ADAFs (e.g., Abramowicz et al. 1988; Chen & 
Taam 1993). The character of the flow in this region is mainly determined by the Euler equation 
(Eq. [2.2]). This equation can be approximately rewritten as: 



(5.1) 



For r 2> r s , the right-hand-side of this equation is always negative since ft < Uk and the third 
term is negligible. If the flow becomes super-Keplerian (in the case of small a) the rhs will become 
positive so that P has a maximum. Even if the flow remains sub-Keplerian (in the case of large 
a), our calculations show that the rhs will still become positive since the third term becomes 
important near r = r s where (H/r) 2 m 1 and M > 1. Therefore P has always a maximum. The 
equation for pressure p is 



1 dp 

p dr 



(JL) 2 _ 1 + ( ^_ 1)7 ( 



(5.2) 



so that for globally sub-Keplerian flows dp/dr is always negative and a pressure maximum is 
present only when the flow becomes somewhere super-Keplerian. 

The topology of the solutions describing rotating accretion flows on to a black hole (Chen et 
al. 1995) implies a different behaviour of the optically thin and optically thick ADAFs at large 
radii. 
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In the optically thick case, the thermal equilibria of the disk at a given radius form a 
characteristic S-shaped curve on the £ — rh plane. Optically thick ADAFs (i.e., slim disks) 
correspond to the upper branch of the S-curve (see Abramowicz et al. 1988; Chen & Taam 1993). 
For a fixed accretion rate the S-curves move up (towards higher rh) with increasing radius (e.g., 
Fig. 2. of Chen & Taam 1993), so that at sufficiently large radii the rh = const solution will 
correspond to an optically thick, gas-pressure dominated, Shakura-Sunyaev solution. Between the 
inner ADAF and the outer stable Keplerian solution there will be, however, an unstable region 
corresponding to the middle branch of the S-curve, so that such global solutions are thermally and 
viscously unstable. 

In the optically thin case the solution structure is different. Contrary to the optically thick 
case where each value of the accretion rate corresponds to exactly one solution, the optically thin 
disk has two solutions for one given rh < m max : an optically thin ADAF and a locally cooled 
optically thin solution. 

Strictly speaking, the self-similar solution can never be an asypmtotic solution for an ADAF 
(at least in the case of bremsstrahlung cooling) because in a self-similar solution Q + = Q a dv oc r~ 3 
while Q- oc r~ 5 / 2 . For large values of a the relevant radius at which local cooling becomes 
important is too large to be of interest but for small viscosity parameters this radius could be 
r < 10 3 — 10 4 as it was shown in the previous section. 

The dimensionless advection factor, £, is shown to be dependent on the fraction of the 
advective cooling, /, since for / < 1 the solution is no longer self-similar. In global solutions, 
the relation between £ and / is nonlinear, but approximately, the smaller the / the smaller the 
£. Furthermore, the value of £ calculated from the self-similar solution is also not valid in the 
transonic region even if / ~ 1. This is because a transonic solution is not self-similar. 

In general, for large enough radii, the solution ceases to exist due to the accretion rate of the 
model exceeding the local maximum rate allowed for such kind of flows. 

Solutions describing optically thin ADAFs can always be, formally, asymptotically Keplerian. 
In Figure 4 we show solutions in which the outer radius was fixed at either r = 2500 r g (heavy 
lines) or r = 9900 r g (thin lines), and the disk parameters are a = 0.01 and rh = 10~ 5 (Model 3). 
The solid lines represent the solution in which a self-similar solution (totally advection-dominated) 
is imposed at the outer boundary; and the dashed lines represent the solution in which a Keplerian 
local type advection-dominated solution is imposed at the outer boundary. Since ADAFs are 
asymptotically self-similar, the Keplerian boundary condition is rather far from the angular 
momentum corresponding to a self-similar flow. That is why the asymtotically Keplerian solution 
takes longer to relax than the asymtotically self-similar one which relaxes in a thin transition zone. 



6. CONCLUSIONS 
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Our calculations of global transonic solutions describing ADAFs show that the self-similar 
solution is generally a good approximation of the full solution as long as the flow is fully 
advection-dominated and far away from the transonic region. Only when the accretion rate is 
close to the maximum one, the solution, as expected, ceases to be advection-dominated and a 
simple self-similar solution is nowhere close to the full solution. 

In some applications (e.g., Narayan et al. 1996; Lasota et al. 1996b) one requires the presence 
of a cold Keplerian accretion disk at some distance r > 10 — W A r g . The transition radius between 
the two types of flows will be given by physical processes that lead to the formation of an ADAFs 
(Meyer & Meyer-Hofmeister 1994; Narayan et al. 1996) This problem was not considered in the 
present work but the results presented above are a step towards its solution. 

Ramesh Narayan, Shoji Kato and Fumio Honma obtained independently global solutions for fully 
advection dominated accretion flows. We thank Andrew King for inspiring remarks. 
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Fig. 1. — Disk solutions for Model 1 (a = 0.1 and m = 10~ 5 , solid lines) and Model 2 (a = 0.1 
and m = 10~ 2 , dashed lines). The heavy lines are the global solution and the thin lines are the 
corresponding self-similar solution. Note that M, f2, and £ depend on m very weakly. 



Fig. 2.— Disk solutions for Model 3 (a = 0.01 and m = 10 -5 , solid lines) and Model 4 (a = 0.001 
and m = 10~ 5 , dashed lines). The heavy lines are the global solution and the thin lines are the 
corresponding self-similar solution. Note the super-Keplerian angular momentum and the maximum 
pressure near the transonic region. Note also that in Model 4, the local cooling becomes important 
for large radii. 



Fig. 3. — Disk solutions for Model 4 (a = 1.0 and m = 10~ 3 ). The heavy lines are the global 
solution and the thin lines are the corresponding self-similar solution. This solution is very sub- 
Keplerian and the sonic point is a nodal type point. 



Fig. 4. — The outer boundary effects shown for Model 3 (a = 0.01 and rh = 10~ 5 ). The heavy and 
thin lines are for the outer radius fixed at r = 2500 r g and r = 9900 r g respectively. The solid lines 
represent the solution in which a self-similar solution (totally advection-dominated) is imposed at 
the outer boundary; and the dashed lines represent the solution in which a Keplerian local type 
advection-dominated solution is imposed at the outer boundary. 



-13- 



Model a M / Me i-in/^K^fg) r s /r g sonic point 

1 0.1 10" 5 0.92967 2.47 saddle 

2 0.1 10" 2 0.92983 2.47 saddle 

3 0.01 10" 5 1.0500601 2.112 saddle 

4 0.001 10" 5 1.0740025 2.048 saddle 

5 1.0 10~ 3 0.446 7.82 nodal 



Table 1: Model Sequences 



